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Given  herein  is  an  easily  implemented  method  for  obtaining  from  complete 
or  censored  data,  approximate  tolerance  intervals  associated  with  the  uppe*  tail 
of  a Weibull  distribution.  These  approximate  intervals  are  based  on  point 
estimators  that  make  essentially  most  efficient  use  of  sample  data.  They  agree 
extremely  well  with  exact  intervals  (obtained  by  Monte  Carlo  simulation  proce- 
dures) for  sample  sizes  of  about  10  cr  larger  when  specified  survival  proportions 
are  sufficiently  small.  Ranges  over  which  the  error  in  the  approximation  is 
within  2 percent  are  determined. 

The  motivation  for  investigation  of  the  methodology  for  obtaining  the 
approximate  tolerance  intervals  was  provided  by  the  new  formulation  of  Lanchester 
Combat  Theory  by  Grubbs  and  Shuford,^  which  suggests  a Weibull  assumption  for 
time-to-incapacitation  of  key  targets.  With  the  procedures  investigated  herein, 
one  can  use  (censored)  data  from  battle  simulations  to  obtain  confidence  intervals 
on  battle  times  associated  with  given  low  survivor  proportions  of  key  targets 
belonging  to  either  specified  side  in  a future  battle.  It  is  also  possible  to 
calculate  confidence  intervals  on  a survival  proportion  of  key  targets  correspond- 
ing to  a given  battle  duration  time. 


This  research  was  supported  by  the  Army  and  the  Navy  through  the  Office  of  Naval 
Research  under  Contract  N00014-73-C-0474  (NR042-321). 
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I.  INTRODUCTION 

In  the  more  usual  analyses  of  life  data,  the  parameters  typically  of 
greatest  interest  are  distribution  percentiles  in  the  lower  tail  of  the  asso- 
ciated life  distributions.  These  percentiles,  in  reliability  analysis,  are 
values  of  reliable  life  corresponding  to  specified  high  survival  probabilities. 

In  contrast,  if  a probability  distribution  is  assumed  for  time-to-inaapaoitation 

of  key  targets  as  in  the  new  formulation  of  Lanchester  Combat  Theory  by  Grubbs 
f3l 

and  Shuford,.  then  it  is  very  often  the  upper-tail  percentiles  of  the  distri- 
bution that  are  of  primary  concern.  There  are,  doubtless,  other  life-time 
distributions  in  which  times  associated  with  low  survival  probabilities  are 
often  of  concern.  For  example,  this  might  be  true  of  distributions  of  particular  v 
types  of  medical  data  or  distributions  of  data  collected  to  demonstrate  inferior 

* 

quality  in  a material  or  piece  cf  hardware. 

Here  we  assume  a two-parameter  Wei  bull  model  for  time-to-failure  or  time- 
to- incapacitation,  as  suggested  by  Grubbs  and  Shuford, ^ and  we  consider  the 
problem  of  obtaining  confidence  intervals  for  distribution  percentiles  asso- 
ciated with  the  upper  tail  of  the  distribution,  or  specified  low  survival  pro- 
portions. Alternatively  a "mission  time"  might  be  specified  and  a confidence 
interval  desired  for  a corresponding  (low)  survival  proportion.  For  Landtester 
Combat  Theory  such  confidence  intervals  could  be  obtained  from  oata  generated 
during  computer  simulation  of  a sample  engagement  that  is,  as  Grubbs  and  Shuford^ 
remark,  "representative  of  the  hypothesized  general  characteristics  of  many 
battles  in  the  supposed  environment." 

"In  particular,  for  example,  we  may  be  interested  in  running  a sample  simu- 
lation of  a combat  situation  in  order  to  see  whether  or  not  it  is  likely  that 
our  choice  of  weapons,  the  tactics  employed  in  using  them,  and  certain  conmand- 
and-control  principles  would  overwhelm  and  defeat  an  eneny  with  somewhat  differ- 
ent weapons  capabilities  in  the  same  hypothesized  battle  environment." 
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If  a Weibull  model  ts  assumed  for  ti.me-to-incapacitation,  or  time- 
to-neutrallze,  for  key  targets  of  each  of  the  two  opposing  sides  in  a fixed 
combat  situation,  then  the  simulated  "observed"  times  of  incapacitation  of 
targets  on  each  of  the  two  sides  can  be  used  to  estimate  the  two  sets  of 
distribution  parameters-  If  a survival  proportion  is  specified  for  either 
side,  then  a time  corresponding  to  this  proportion  can  be  estimated  as  a func- 
tion of  the  parameter  estimates.  If,  alternatively,  a "mission  time,"  or 
"battle  duration  time,"  is  specified,  corresponding  survival  proportions  can 
be  estimated  for  both  sides  from  the  two  sets  of  parameter  estimates. 

In  the  following,  an  approximation  that  can  be  used  in  conjunction  with 
the  parameter  estimates  for  obtaining  confidence  intervals  on  survival  pro- 
portions corresponding  to  specified  battle  duration  times  or  on  duration  times 
corresponding  to  specified  survival  probabilities  is  described. 

II.  A WEIBULL  MODEL  FOP  TIME  TO  INCAPACITATION 

It  is  assured  here  that  the  random  variate  time-to-iru^cu?Ltation  T is  such 

that 


(1  - exp[-(t/6)6]  , t > 0 

P(T  < t)  * ! 


I 


(1) 


otherwise. 


where  6,  8 > 0.  For  X 3 fcnT,  n = 2-n6  and  £ 3 S ^ > 0,  this  is  equivalent  to 


P(X  < x)  3 1 - expj-exp[(x-n)/£]j 


(2) 


The  random  variati  T has  a Wei  bull  distribution  with  shape  parameter  8 and  scale 
parameter  ( characteristic  time-to-inccnpaeitation ) 6.  (The  parameters  8 and  6 
will  of  course  tend  to  be  different  for  the  two  opposing  sides  in  a battle.) 
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The  random  variate  X has  the  first  asymptotic  distribution  of  the  smallest 
extreme  (the  extreme-value  distribution).  The  parameter  n is  the  mode  and 
ir§//6  is  the  standard  deviation  of  the  distribution  of  X.  The  linear  estimates 
discussed  in  the  following  are  based  on  observations  on  X ordered  from  smallest 
to  largest. 

C31 

Grubbs  and  ShufordL  J reference  three  papers  containing  tabulations  from 
which  one  can  calculate  confidence  intervals,  for  each  of  the  two  opposing 
sides,  on  the  distribution  parameters,  on  high  survival  probabilities  of  key 
targets,  or  on  lower-tail  distribution  percentiles.  These  are  based  on  (itera- 
tively obtained)  maximum-likelihood  estimates,  best  linear  invariant  (BLI) 
estimates,  or  approximations  of  BLI  estimates.  One  of  these  three  papers 
(Bill man,  Antle  and  Bain‘-,J)  applies  only  to  censored  samples.  The  tables  given 
are  extensions  of  those  published  bv  Thoman,  Bain  and  Antle^^  applying  to 
maximum-likelihood  estimates  obtained  from  complete  samples.  The  paper 
referenced  in  Grubbs  and  Shuford  with  tables  or  obtaining  Weibull  confidence 
intervals  from  best  linear  invariant  estimates  of  MannL'  * J is  by  ann  and 
Fertig.^1^  The  tabulations  given  therein  apply  to  samples  of  sizes  2 through 

25,  either  comp! ete  or  with  all  possible  right-hand  censorings.  Tabulations 

m 

of  Johns  and  LiebermanL  J pertain  to  asymptotic  approximations  of  best  linear 
invariant  estimates  and  four  right-hand  censorings  for  each  of  six  sample  sizes 
ranging  from  10  to  100.  These  published  tabulations  have  all  been  generated  by 
means  of  Monte  Carlo  simulation  techniques  since  the  exact  distributions  of  the 
estimators  cannot  be  determined.  The  estimators  to  which  they  apply  are  all 
asymptotically  fully  efficient;  that  is,  as  sanple  si  :e  n approaches  infinity 
the  variances  of  the  (asymptotically  unbiased)  estimators  approach  the  Cramer-Rao 
lower  bounds  for  variances  of  regular  unbiased  estimators. 


I 1 

j 

Ij 
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The  parameter  of  primary  int.=  est  in  the  present  context  is  tp  * exp(xp), 

the  10QP th  percentile  of  the  distr  bution,  where  P is  proportion  of  key  targets 

* 

incapacitated  by  time  tp.  The  proportion  P incapacitated  will,  of  course,  tend 
to  be  different  for  each  of  the  two  opposing  sides  for  fixed  tp.  Note,  from 
(2),  that  xp  * n + 5n[-&n(l-P)].  If  P is  specified,  then  essentially  optimum 
point  estimates  of  tp  in  terms  of  mean  squared  error  in  the  log  space  can  be 
obtained  for  targets  on  each  of  the  two  sides  as  expjn  + ££n[-£..-:l-P)]}  or 
expjn  + ££n[-£n(l-P)]} , where  for  each  side  n and  £ are  best  linear  invariant 
estimates  and  n and  £ are  maximum- likelihood  estimates  of  n and  £,  respectively, 
and  Zn[-£n(l-P)]  is  the  lOOPth  percentile  of*(X-n)/£.  If  a battle  duration 
time  tQ  * exp(xQ)  is  specified,  then  exp{-exp[(x0  - n)/£]f  or 
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! 
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exp{-exp[(xQ  - n)/£]}  provide  efficient  point  estimates  of  l-Pfty),  the  propor- 
tion of  targets  on  each  side  surviving  incapacitation  until  at  least  time  tQ. 

As  noted,  the  tabulations  mentioned  earlier  apply  only  to  calculation  of 
confidence  bounds  on  lower-tail  distribution  percentiles  or  on  high  survival 
probabilities-  Because  of  this  and  because  these  tabulations  of  exact  values 
were  necessarily  generated  by  Monte  Carlo  simulation  procedures,  one  would 
expect  that  similar  simulation  procedures  would  be  necessary  in  order  to  find 
values  to  be  used  in  calculating  confidence  bounds  on  upper-tail  distribution 
percentiles  or  low  survival  proportions.  We  show  now  that  such  is  not  the  case. 

III.  APPROXIMATIONS  FOR  CONFIDENCE  INTERVALS  ON  BATTLE  OURATIQN  TIME 

Recently,  Engelhardt  and  3ain^  and  Mann  and  Fertig^^  have  shown  that 
efficient  linear  estimators  of  £ (those  with  smallest  or  nearly  smallest  mean 
squared  error)  are  approximately  proportional  to  chi-square  variates.  Thus, 

pi  w*! 

using  the  two-moment  fit  of  Patnaik,“  one  sees  that  if  £*  n is  an  unbiased 

efficient  linear  estimator  of  £,  based  on  the  r smallest  observations  from  a 

sample  of  size  n,  with  variance  C £ , then  2£*  / (C_  „£)  is  approximately  a 

i jii  r ,n  r *n 

chi-square  with  2/Cr  n degrees  of  freedom. 


In* «*»« 


If  i*  „ is  the  unique  best  {uniformly  minimum  variance)  linear  unbiased  estimator  . 

• *n 


of  5,  then  I a £*  _/(l+C_  _)  is  the  best  linear  invariant  estimator  of  £. 

r *«*  r *n  r *n 


(The  estinetor  £„  has  mean  squared  error  independent  of  n and  uniformly  smaller 

• »n 


mean  squared  error  than  that  of  . ) As  can  be  seen  from  the  results  of  Lawless 

r jii 

and  Mann,^  is  very  nearly  equivalent  to  the  iteratively  obtained  maximum- 


likelihood  estimator  % of  £.  This  near  equivalence  can  be  seen,  too, 

r*n 

if  one  compares  tabulations  of  distribution  percentiles  of  £ /£  and  of 

r ,n 


r,n 


appearing  in  Mann  and  Fertig^^  and  Thoman,  Bain  and  Antle,^8^  respectively. 


keeping  in  mind  that  both  sets  of  Monte  Carlo  generated  tabulations  are  correct 
to  within  about  a unit  in  the  second  decimal  place  for  the  values  of  n compared. 
Hence  2(1+C  )\  /( C £)  and  2(1+C  )%  _/(C  £)  are  both  approximate  chi- 

* i*»  * jii  r jii  r fit  i 9 n r fil 

square  variates  with  2/C„  „ degrees  of  freedom.  This  fact  allows  one  to  calcu- 


r,n 


late  confidence  bounds  for  £ by  using  either  iteratively  obtained  maximum- 
likelihood  estimates  or  best  linear  unbiased  or  best  linear  invariant  estimates 
(or  efficient  approximations  thereof)  in  conjunction  with  values  of  C 


r,n 


obtainable  from  tabulations  in,  for  example,  Mann,-7*8-^  Engelhardt  and  Bain,^ 


Mann,  Schafer  and  Singpurwalla,^4^  and  Mann  and  Fertig.^2^  Engelhardt  and 


Bain  also  provide,  for  various  values  of  r/n,  approximate  expressions  for  C 

r,n 

that  are  quadratics  in  1/n. 

Mann,  Schafer  and  Singpurwalla^4^  use  the  chi-square  property  of 
25r,n^Cr,n^  t0  constnjct  an  approximately  F-distributed  statistic  for  obtaining 
confidence  intervals  on  tR  (reliable  life,  or  in  this  context.,  battle  duration 
time).  First,  consider  best  (or  approximately  best)  linear  unbiased  estimators 


iiiiiiiiuui 


Define  Xr,n  as  K,n  ’ (Br,n/Cr,n>?r,JH^r,n  ' (Br.n/Cr.n)gr,J,rith  var1ance 
' Br,n/Cr,n**  The  covariance  of  Xr,n  ard  n is  given  by 
[3r  n - (Br  n/Cr  n)Cr  n]£2  = 0.  Mann,  Schafer  and  Singpurwalla^1^  give  examples 

of  comparisons  with  exact,  Monte  Carlo-generated,  previously-tabulated  values  to 
show  that  for  small  values  of  P one  can  consider 


(Xp,n  " xp^|i*Br,n/,f''r,n  " ^C-^>(l-P)]}Cr>r 


to  have  an  F distribution  with 


- 2|WCr,r.  + “C-in(l-P)]p/(Ar_n  - B2>(,/Cr>n) 


v2  = 2/Cr,n 


♦ 1 
degrees  of  freedom.  Note  that  both  C*  n and  (X*  n - xp)/{-Br  n/C„.  n - in[-£n(l-P)]( 

are  unbiased  estimators  of  5;  and  with  a two-moment  fit  to  chi-souare,  each  is 
under  the  proper  conditions,  an  approximate  chi-square  over  its  degrees  of  free- 
dom. The  value  of  is,  like  the  value  of  Vg,  obtained  from  the  two-moment 

it  a 

chi-square  fit  [2m(X  - xp)/v  is  approximately  chi-square  witn  2m/v  degrees  of 

freedom,  for  E{X*  - xp)  = m and  var(X*  - xp)  = v]  and  can  be  calculated  from 

tabulated  values.  Values  of  A„  _ and  B„  _ can  be  found  with  tabulations  of 

r,n  r,n 

C„  Additional  values  of  A„  for  large  n,  are  given  by  Mann,  et  al,^^ 

r y n i 91I 

P.  252.  X 
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Mann^^  gives  the  rationale  for  the  chi-square  fit  for  the  numerator  arid 
gives  ranges  in  terms  of.  v.  and  Vg  fcr  which  values  for  obtaining  confidence 
bounds  on  Xp  from  the  F-approximation  (3)  are  in  error  by  two  percent  or 
ie;s  for  P * .01,  .05  and  .10.  In  most  cases  the  approximation  is  excellent, 

f §1 

even  for  n or  r as  small  as  4 unless  cen'orlng  is  extreme.  Law! ess c J shows 

that  lower  confidence  bounds  on  x0  based  on  (3)  with  X*  _ replaced  by 

r r *n 

n - (B„  JC  . a function  of  the  maximum-likelihood  estimators  n_  „ , 

r,n  r,n  r,n  r»n  r,n 

and  tr  n»  and  £*  n replaced  by  {1  + Cr  n)€r  n agree  to  two  or  more  significant 

figures  with  exact  values  for  sample  sizes  25  through  60  unless  nearly  90  per- 

rq'i 

cent  of  the  sample  is  censored;  i.e.,  r/n  < .10.  MannL  J demonstrates  the 
excellence  of  a similar  F-approximation  in  obtaining  prediction  intervals  for 
future  samples,  or  lots  to  be  manufactured  in  the  future. 

Thus,  for  a specified  large  survival  proportion  P applying  to  either  of  the 
opposing  sides  in  a-’battle  one  might  conclude  from  (3)  that  an  approximate 
lower  (l-a)-level  confidence  bound  on  the  corresponding  battle  duration  time 
tp  * exp(xp)  can  be  calculated  as 


e*plXr,i»  " Fa(VVHBr,n/Cr,n  + ^"(l-P)]^ 


where  Fj(vi .v^)  is  the  Iddath  percentile  of  F with  v-j  and  v 2 degrees  of  freedom 
and  and  v2  defined  by  (4)  and  (5),  respectively.  An  upper  (l-a)-level  confi- 
dence bound  on  tp  is  given  by 

HC  + Fl-a(VV(E,r,n/Cr,n  + *nC-*»C1-P)3K*}  <?i 


A two-sided  confidence  interval  on  tp  at  confidence  level  (l-2a)  is  an  interval 
with  lower  and  upper  bounds  g*ven  by  (6)  and  (7),  respectively. 
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Recall,  for  use  of  f tables,  that  ,v2)  = l/Fj^Vg,^).  In  general, 
however,  and  Vg  will  not  be  integers  so  that  one  can  interpolate  in  tables 
of  percentiles  of  F,  or  can  use  *he  approximation  suggested  by  Mann  and 

f13l 

Grubbs*  J to  evaluate  F (v, -v?)  for  given  y. 

It  should  be  clear  from  inspection  of  the  expression  (3),  (4)  and  (5)  that 
if  xp  is  specified  and  a confidence  interval  is  required  for  P (or  equivalently, 
for  in[-in(1-?}]),  it  must  be  confuted  iteratively  since  P occurs  *0  both  (3) 
and  (4).  In  this  case  we  find  that  a lower  confidence  bound  for  P at  levs! 
l-o  for  specified  xp  = ln(tp)  is  given  by  1 - expj-exp{(xp  - X*  n)/ 

[5r  nFi-0(y,j»v2^  * Br  r/cr  n}|»  where  is  a function  of  the  value  of  P which 
is  being  determined. 

IV.  PRECISION  OF  THE  APPROXIMATION 

A Monte  Carlo  simulation  study  was  undertaker,  so  that  the  precision  of  the 
F-approximaticn  (3)  could  be  determined  for  values  of  P of  Interest.  The  values 
of  P considered  were  .75(.05}.95  and  .99  and  the  percentiles  of  the  appropriate 
F-distributions  tabulated  and  compared  with  Monte  Carlo  values  were  .01,  .02, 

. 05 ( . 05 ) . 95 . .93,  .39.  The  values  of  n and  r considered  were  $!  * 8,  r * 4,8 
and  n * 15,  r * 5,10,15.  For  r/n  fixed,  one  would  expect  precision  to  increase 
with  increasing  sample  size  r;  since  the  basis  of  the  approximation  is  an 
asymptotic  result.  Monte  Carlo  sample  size  for  the  study  was  10,000. 

It  was  found  that  the  precision  of  the  approximation  increased  with 
increasing  P over  the  range  considered  and  with  a limited  amount  of  increased 
censoring  for  a fixed  sarole  size.  In  both  cases  the  increase  in  precision 
resulted  apparently  from  an  increase  in  v,  relative  to  v„.  Results  indicate  that 

* C 

agreement  of  exact  percentiles  of  (3)  and  those  based  on  the  F-approxiroation 
is  to  within  about  2 percent  over  the  range  of  percentiles  f**om  .01  to  .99  for 
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mipit# 


v-j  > .03v2  + 20.0  , > 8.0 


and  over  the  range  of  percentiles  from  .10  to  .99  for 
v1  > 0.3v2  + 4.0  , V2  > 8.0  . 

These  ranges  are  similar  in  spirit  to  those  given  by  Mann~10^  for  small  values 
Ot  P.  Shown  in  Table  I are  the  results  of  six  of  the  3C  independent  simulations 
performed 

V.  EXAMPLE 

An  example  of  simulated-battle  data  is  giver,  in  Table  I of  Grubbs  and 
Shuford,^  and  a description  of  the  simulated  battle  is  given  preceding  their 
table.  For  each  of  the  two  sides  in  the  engagement,  the  number  of  key  targets 
(tanks,  in  this  case)  is  20.  The  data  consist  of  times-to-incapacitation  in 
minutes  for  four  CBT's  (chief  battle  tanks)  and  for  five  RIO  tanks  on  the 
opposing  side.  Extensions  of  tables  of  Mann^  have  been  used  to  provide  esti- 
mates n4  go  * 5,327  and  20  = ^-002  from  the  RIO  data. 


To  use  (7)  to  obtain  a confidence  interval  on  tp  for  the  RIO  tanks  for  a 
specified  P of,  say,  .95,  one  needs  to  calculate  in[-in(l-P)]  * 1.0972  and  to 
look  up  tabulated  values  from  which  values  of  A,  on,  3-  on  ar.d  C-  can  be 

3,ZU  b,£U  3,£U 

determined  in  the  report  that  provides  the  coefficients  (or  weights) 

for  calculating  the  estimates  of  n and  £,  The  tabulated  values  from  which  the 

necessary  constants  can  be  calculated  are  Eg  20'^  3 £5  20^n~n' 

E5,20(LB)  * £5,2oC^-5)2]/52  and  E520(CP)  * £C(n-n)(C-C)3/?2.  To  determine 
Ag  2g»  B-  2q  and  ^5  20  one  nse<*s  fcncw  relationships: 


lllllllll^llllllll^Diidlim " 1MMA1AL1! sSMaritii M ilslta nWltow M ’ Vt  i!H.  0*1 


Table  I.  Approxinate  and  Monte  Carlo  (M.C.)  Values  of  lOOrt h Distribution 
Pe**centiles  of  the  Approximate  F-Variate  (3) 


n * 15,  r * 10,  P * .85 
vj  51  15.6  v2  = 22.4 


n » 15*  r - 5,  P * .85 
vt  = 23.2  v2  = 3.8 


r.  - 15,  r * 5,  P 
Vj  = 19.1  v2  * 


M.C. 

Percentile 


Approximate 

Percentile 


0.268 

0.318 

0.406 

0.491 

0.697- 

0.891 

1.04 

1.21 

1.55 

2.30 

2.96 

4.18 

5.14 


0.329 

0.412 


0.707 


n = 15,  r = 

15,  P = .99 

vi  * 40.4 

v2  = 44.1 

M.C. 

Approximate 

T 

Percentile 

Percentile 

0.01 

0.472 

0.482 

0.02 

0.05 

0.510 

0.588 

0.10 

0.665 

0.670 

0.25 

0.808 

G.809 

0.40 

0.922 

0.922 

0.50 

0.996 

0.998 

0.60 

1.07 

1.08 

0.75 

1.23 

1.23 

0.90 

1.49 

1.49 

0.95 

1.68 

1.67 

0.98 

1.93 

1.90 

0.99 

2.10 

2.07 

n * 8,  r * 4,  P 
vj  * 23.0  v2 


M.C. 

Percentile 


0.284 

0.322 

0.407 

0.498 

0.700 

C.914 

1.07 

1.26 

1.71 

2.63 

3.55 

5.43 

7.20 


n = 8,  r * 
vi  = 9.9 

8.  P * .95 
v2  * 21.5 

M.C. 

Percentile 

Approximate 

Percentile 

0.059 

0.147 

0.297 


0.652 


0.971 

1.12 

1.40 


and 


A5,20  = E5,20^  + ^E5,20^CP^  ^ " E5,20^LB^ 

From  the  tables  and  these  relationships,  we  find  Ag  2g  = 0*70308,  Bg  2g  = 0.33548 

and  Cg  20  = 0.23662.  (For  n = 25(5)'  \ p = r/n  = 0.10(0.10)1.00,  values  of 

A , B and  C (or  l „)  are  tabulated  directly  by  Mann,  Schafer  and 
r,n  r,n  r,n  r,n 

Singpurwalla,^4^  p.  252  and  p.  244.) 

Next,  from  (4)  and  (5),  we  find  for  the  RIO  data. 


v1  = 2(1.4178  t 1 .0972) Z/ (0.2274)  = 55.62 


v„  = 2/. 23662  = 8.45  , 

C. 


correct  to  two  decimal  places.  Note  that  v2  > 8.0  and  v-j  > .3(v2)  + 20  * 22.53.  * 

so  we  might  expect  the  F-approximation  to  give  very  nearly  the  correct  percentile 
of  the  distribution  of  (3)  over  the  percentile  range  of  .01  through  .99.  Using 
the  method  of  Mann  and  Grubbs  to  evaluate,  say,  the  7 0th  percentile  of 
F(55.62,8.45),  we  obtain  0.567.  Then  a 90  percent  lower  confidence  bound  on 
t g5  is,  from  (6),  given  by 

exp[5.040  - 1.4178(1.002)  + 0.567(1.4178  + 1. 09721(1. 23662)(1. 002)]  = 218  minutes 


If  a lower  confidence  bound  on  6 is  desired,  as  on  page  938  of  Grubbs  and 
Shuford,^  ther.  the  specified  value  of  P is  1 - exp(-l)  = .63  and  £n[-in(l-P)]  * 0 
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tiiim 


so  that  v.]  a 2(1.4178)  v(.2274)  = 17.68  and  Vg  = 8.45,  as  before.  The  value 
of  F 05(17.68,8.45)  is  approximately  0.408  so  that  a 95  percent  lower  confidence 
bound  on  6 = exp(n)  is  therefore  given  by 

exp [5. 040  - 1.4178(1.002)  t .408(1 .4178)(1 .23662)0 .002)]  = 76.4  minutes  . 


This  agrees  welt  with  the  lower  95  percent  confidence  bound  74.4  minutes  shown 
in  the  example  of  Grubbs  and  Shuford  and  calculated  as  exp(n  - w ggC).  where 
w gj.  is  a percentile  of  the  distribution  of  (n-n)/5  computed  by  Monte  Carlo 
simulation  procedures  and  published  in  the  report  by  Mann,  Fertig  and  Scheuer 
with  tabulations  supplemental  to  those  of  Mann  and  Fertig. A lower  confi- 
dence bound  on  6,  incidentally,  will  correspond  to  a lower  confidence  bound  on 
mean-tiMe-to-inccrpacitation  for  8 = £ - 1 * since  the  exponential  distribution 
is  a special  case  of  the  Wei  bull  with  shape  parameter  equal  to  1. 
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